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Abstract 



We have recently developed a coupled chemistry-emission 
model for the green (5577 A) and red-doublet (6300, 6364 
OA) emissions of atomic oxygen on comet C/1996 B2 Hyaku- 
take. In the present work we applied our model to comet 
__! C/1995 01 Hale-Bopp, which had an order of magnitude 
higher H2O production rate than comet Hyakutake, to eval- 
uate the photochemistry associated with the production and 
'Z. loss of 0( 1 S) and 0( 1 D) atoms and emission processes of green 
r-joand red-doublet lines. We present the wavelength-dependent 
photo-attenuation rates for different photodissociation pro- 
*£h cesses forming 0( 1 S) and 0( 1 D). The calculated radiative effi- 
'"pTbiency profiles of 0( 1 S) and 0( 1 D) atoms show that in comet 
O Hale-Bopp the green and red-doublet emissions are emitted 
^mostly above radial distances of 10 3 and 10 4 km, respectively. 
^The model calculated [OI] 6300 A emission surface bright- 
i_ness and average intensity over the Fabry-Perot spectrome- 
ter field of view are consistent with the observation of Mor- 
r ^genthaler et al. (2001), while the intensity ratio of green to 
qq red-doublet emission is in agreement with the observation of 
Zhang et al. (2001). In comet Hale-Bopp, for cometocentric 
C> distances less than 10 5 km, the intensity of [OI] 6300 A line is 
*0 mainly governed by photodissociation of H2O. Beyond 10 5 km, 
rHOf'D) production is dominated by photodissociation of the 
water photochemical daughter product OH. Whereas the [OI] 
^1 5577 A emission line is controlled by photodissociation of both 
t 7"!h 2 and CO2. The calculated mean excess energy in various 
^photodissociation processes show that the photodissociation 
of CO2 can produce 0( 1 S) atoms with higher excess velocity 
_j compared to the photodissociation of H 2 0. Thus, our model 
calculations suggest that involvement of multiple sources in 
the formation of 0( 1 S) could be a reason for the larger width 
of green line than that of red-doublet emission lines observed 
in several comets. 



1 Introduction 

Owing to its very high H 2 production rate, C/1995 01 Hale- 
Bopp was a great comet in the night sky of the year 1997. 
The visible emissions of atomic oxygen ([OI] 6300, 6364, and 
5577 A), which are accessible to the ground-based optical in- 
struments, have been observed by Morgenthalcr et al. (2001) 
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and Zhang et al. (2001) in the coma of Hale-Bopp. Since the 
lifetime of oxygen atom in the 1 D metastable state is rela- 
tively small (^110 s) compared to the photochemical lifetime 
of H 2 (~8 x 10 4 s), it cannot travel larger distances in the 
coma without decaying to the ground 3 P state. Moreover, 
most of the production of oxygen in the *D state is through 
photodissociative excitation of H 2 (Bhardwaj and Haider, 
2002). Thus, [OI] 6300 A emission has been used to trace 
the spatial distribution as well as to quantify the production 
rate of H 2 in several comets (Delsemme and Combi, 1976; 
Delsemme and Combi, 1979; Fink and Johnson, 1984; Schultz 
et al., 1992; Morgenthaler et al., 2001; Furusho et al., 2006; 
Fink, 2009). 

Based on the study of Fcstou and Fcldman (1981) the inten- 
sity ratio of green (5577 A) to red-doublet (6300 A and 6364 A) 
emissions (here after G/R ratio) of atomic oxygen has been 
used to determine whether the parent source of these lines 
is H 2 or C0 2 /CO in the coma of comets (Cochran, 1984, 
2008; Morrison et al, 1997; Zhang et al., 2001; Cochran and 
Cochran, 2001; Furusho et al., 2006; Capria et al., 2005, 2008, 
2010; McKay et al., 2012a,b). The modelling studies of these 
emissions in comets showed that photodissociative excitation 
of H 2 is the major production process of the [OI] 6300 A 
emission (Festou and Fcldman, 1981; Bhardwaj and Haider, 
2002; Capria et al., 2005, 2008; Bhardwaj and Raghuram, 
2012). Our recent theoretical study (Bhardwaj and Raghu- 
ram, 2012) for these prompt emissions of atomic oxygen in 
comet C/1996 B2 Hyakutake showed that more than 90% 
of the 0( 1 D) is populated via photodissociative excitation of 
H 2 and the rest through photodissociation of other oxygen 
bearing species, like C0 2 and CO. It also showed that quench- 
ing by H 2 is the major loss mechanism of 0( X D) up to come- 
tocentric distances of 1000 km, and above that distance radia- 
tive decay takes over. The study of Bhardwaj and Raghuram 
(2012) demonstrated that the G/R ratio depends not only on 
the photochemistry involved in populating 0( 1 S) and 0( 1 D) 
atoms in the cometary coma, but also on the projected area 
observed for the comet, which is a function of slit dimension 
used for observation and geocentric distance of the comet. 

In the present study we applied our coupled chemistry- 
emission model (Bhardwaj and Raghuram, 2012) to comet 
C/1995 01 Hale-Bopp, which had an order of magnitude 
higher H 2 production rate compared to that of comet Hyaku- 
take, to evaluate the production and loss mechanisms of 0(*D) 
and 0( X S) and generation of green and red-doublet emissions. 
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Our aim in this paper is to study the photo-attenuation in 
comets having high H2O production rates and its implications 
on the photochemistry of metastable 0( 1 S) and 0( 1 D) atoms. 
We compared our model calculated [01] 6300 A emission sur- 
face brightness profile with the observation of Morgenthaler 
et al. (2001). We have shown that the photodissociation of 
H2O mainly controls the formation of 0( 1 D) and subsequently 
determines the intensity of [OI] 6300 A emission. However, 
in the case of [OI] 5577 A emission, the photodissociation of 
both H 2 and CO2 plays an important role in the formation 
of metastable 0( 1 S), with photodissociation of C0 2 being the 
major production source of 0( 1 S) in the inner cometary coma. 
We suggest that in comets with significant (>5%) CO2 rela- 
tive abundance, the photodissociation of C0 2 is more effective 
in populating 0( 1 S) than the photodissociation of H2O. The 
model calculated G/R ratio is consistent with the observed 
value of Zhang et al. (2001). Based on the model results, we 
suggest that the involvement of multiple species in the for- 
mation 0( 1 S) could be a reason for the width of the green 
line emission to be larger than the red-doublet emission lines 
observed in several comets by Cochran (2008). 

2 Model 

The details of the model and the chemistry are presented in 
our previous work (Bhardwaj and Raghuram, 2012). Here 
we present the input parameters that have been used in the 
model for the observed condition of comet Hale-Bopp on 10 
March 1997 (solar radio flux F10.7 = 74.7 x 10~ 22 J s" 1 m~ 2 
Hz -1 ) when it was at a geocentric distance (A) of 1.32 AU 
and a heliocentric distance (r^) of 0.93 AU. The photochemi- 
cal reaction network and cross sections for photon and electron 
impact processes are same as used in the previous work and 
any changes made are mentioned. The degradation of solar 
EUV-generated photoelectrons is accounted by using Analyt- 
ical Yield Spectrum (AYS) technique which is based on the 
Monte-Carlo method (Singhal and Bhardwaj, 1991; Bhardwaj 
and Singhal, 1993; Bhardwaj and Michael, 1999a, b; Bhard- 
waj and Jain, 2009). Details of the AYS approach and the 
method to calculate photoclcctron flux and excitation rates 
are given in our earlier papers (Bhardwaj et al., 1990, 1996; 
Bhardwaj, 1999, 2003; Haider and Bhardwaj, 2005; Bhardwaj 
and Raghuram, 2011; Raghuram and Bhardwaj, 2012; Bhard- 
waj and Jain, 2012). 

In our previous work (Bhardwaj and Raghuram, 2012) it has 
been shown that the contribution of several processes to the 
production of 0( X S) and 0( 1 D) is small. Thus, only important 
production and destruction mechanisms of metastable 0( 1 S) 
and 0( 1 D) are presented in Table 1. The branching ratio of 
0.81 is used for the production of OH in the photodissociation 
of H 2 (cf. Hucbner et al., 1992; Nee and Lee, 1984). The ra- 
dius of the cometary nucleus is taken as 25 km (Weaver et al., 
1997; Combi et al., 1999). Though a variety of O-bearing mi- 
nor species (with relative abundances <1%, Bockelee-Morvan 
et al., 2000) have been found in comet Hale-Bopp, the dom- 
inant species H2O, CO2, and CO are only considered in our 
model calculations. The neutral density profiles of these par- 
ent species are calculated using Haser's formula. 



Using ground-based observations of direct H2O infrared 
emissions during January to May 1997, Dello Russo et al. 
(2000) derived water production rates at different heliocen- 
tric distances and fitted as a function of as Qh 2 o — 8.35 ± 
0.13 x 10 30 [ r ^ 1 - 88±013 )] molecules s" 1 . We used this expres- 
sion as a standard input in calculating H 2 density profiles on 
different days. The H2O production rates in Hale-Bopp are 
also derived by observing the emissions from its dissociative 
products, like OH 18-cm emission, OH (0-0) 3080 A emission, 
[OI] 6300 A emission, and H Lyman-a, over large projected 
distances (Weaver et al., 1997; Colom ct al., 1999; Schleicher 
et al., 1997; Combi et al., 2000; Woods ct al., 2000; Morgen- 
thaler et al., 2001; Harris ct al., 2002; Fink, 2009). These 
derived H2O production rates depend on the observational 
condition and also on the assumed model parameters. We 
have discussed the effect of H 2 production rate on the cal- 
culated green and red-doublet emission intensities of atomic 
oxygen in the Section 4.2.2. 

High resolution ground-based infrared spectroscopic obser- 
vations are used to study the CO production rate in comet 
Hale-Bopp from June 1996 to September 1997 (DiSanti et al., 
2001). The spatial distribution of CO in the coma of Hale- 
Bopp is assumed to have two distinct sources: nucleus- 
originated, and extensively distributed in the cometary coma. 
During this observation period, the relative abundance of CO 
is found to be -25% of H 2 0. DiSanti et al. (2001) fitted the 
observed CO production rate as a function of heliocentric dis- 
tance near the perihelion as Q C o = 2.07 x 10 30 r - 1 - 66±0 - 22 
molecules s _1 , and above 1.5 AU as Qco = 1-06 x 10 30 
r -i.76±o.26 mo i ecu i es s -i since observations of [OI] 6300 
A emission are done when comet was at around 1 AU, we 
have used the former formulation to calculate the neutral CO 
density in our model. DiSanti ct al. (2001) suggested that 
the increase in CO production rate below 1.5 AU is due to 
distributed sources prevailing in the cometary coma. Recent 
study of Bockelee-Morvan et al. (2010) showed that the in- 
frared CO (1-0) rotational vibrational emission lines are opti- 
cal thick in the cometary coma of Hale-Bopp. Based on the 
modelling studies of these emission lines they rejected the idea 
of extended source distribution of CO in comet Hale-Bopp. 
However, our model calculations show that the role of CO in 
determining green and red-doublet emission intensities is very 
small compared to other species, and hence the impact of dis- 
tributed CO source is insignificant on these forbidden emission 
lines. 

The CO2 has been detected in Hale-Bopp by Crovisier et al. 
(1997) in April 1996, when the comet was at heliocentric dis- 
tance of 2.9 AU. Based on the infrared emissions between 2.5 
to 5 /xm, the derived CO2 production rate at 2.9 AU was 1.3 x 
10 28 molecules s -1 , which corresponds to a relative abundance 
of —20% of H2O. Assuming that the photodissociative exci- 
tation is the main production mechanism in populating the 
CO(a 3 n) metastable state, the observed CO Cameron band 
(a 3 n — > X 1 E + ) emission intensity has been used to estimate 
the abundance of CO2 in this comet by Weaver et al. (1997). 
The estimated CO2 abundance is more than 10% when the 
comet was beyond 2.7 AU. However, our model calculations on 
comets 103P/Hartley 2 (Bhardwaj and Raghuram, 2011) and 
lP/Halley (Raghuram and Bhardwaj, 2012) have shown that 
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Table 1: Major production and destruction processes of the 0( 1 S) and 0( 1 D). Photorates are calculated using solar flux on 10 
Apr 1997 (solar minimum period : solar radio flux F10.7 = 74.7 x 10~ 22 J s _1 m~ 2 Hz -1 .) and scaled to 0.92 AU heliocentric 
distance. 
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0.134 


Slanger et al. (2006) 
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1.26 


Wiese et al. (1996) 
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This work 


OH + 


7.01 x 10" 6 


Huebner et al. (1992) 1 


C0 2 4 
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6.2 x 10 -7 


This work 


CO + 


6.0 x 10" 8 


This work 


H 2 0+ 


4.3 x 10 -7 (300/T € ) 5 x 0.045 f 


Rosen et al. (2000) 


CO+ - 


1- e th -> O^D) + others 


5.0 x 10" 8 x (300/T e ) 46 


Mitchell (1990) 


0( J D) 


+ H 2 -> 2 OH 


2.1 x 10 -10 


Atkinson et al. (1997) 


oCd) 


— > 0( 3 P) + hi/ 

\, 6300/4 


6.44 x 10" 3 


Storey and Zeippen (2000) 


0( 1 D) 


— ► 0( 3 P) + hi/ 

v ' KSMA 


2.15 x 10" 3 


Storey and Zeippen (2000) 



^0.045 is the assumed branching ratio for the formation of 0( 1 S) and 0( 1 D) via dissociative recombination of H20 + ion (see Bhardwaj 
and Raghuram, 2012). *Huebner et al. (1992) calculated this rate using theoretical OH absorption cross section of van Dishoeck and 
Dalgarno (1984). ^Huebner et al. (1992) calculated this rate based on experimentally determined OH absorption cross section of Nee 
and Lee (1984). hv : solar photon; eth '■ thermal electron; T e : electron temperature. 



photoelectron impact excitation is the main production mech- 
anism of CO Cameron band emission and not the photodisso- 
ciation of CO2. Assuming that the CO2/CO abundance ratio 
did not vary with heliocentric distance in this comet, Bockelce- 
Morvan et al. (2004) suggested 6% relative abundance of CO2 
when the comet was at 1 AU. We have taken 6% CO2 rela- 
tive abundance with respect to H 2 in the model. However, 
we discuss the impact of CO2 abundance by varying its rela- 
tive abundance on the calculated intensities of green and red- 
doublet emissions. The OH neutral density profile in comet 
Hale-Bopp is calculated by fitting Harris et al. (2002) observed 
OH (0-0) 3080 A resonant scattering emission along the pro- 
jected distance with the Haser's two step formulation. The 
photodissociative excitation rates of OH producing 0( 1 S) and 
0( 1 D) are taken from Huebner et al. (1992) which were de- 
termined using theoretical (van Dishoeck and Dalgarno, 1984) 
and experimental (Nee and Lee, 1984) photoabsorption cross 
sections, respectively. 

There is a clear evidence that in comet Hale-Bopp the ex- 
pansion velocity of neutrals increases with increasing cometo- 
centric distance (Colom et al., 1999; Bivcr et al., 1997; Harris 
et al., 2002). The sources involved in accelerating the neu- 
tral species across the cometary coma is discussed in several 
works (Colom et al., 1999; Combi et al., 1999; Harris et al., 
2002; Combi, 2002). To incorporate the acceleration of these 
neutrals in our model we have taken the velocity profile cal- 
culated by Combi et al. (1999) at 1 AU and used as a input 
in the Haser's density distribution. We also verified the effect 
of expansion velocity on the calculated intensity of green and 
red-doublet emissions by varying its static value between 0.7 
to 2.2 km s _1 , which is discussed in Section 4.2.1. 



The input solar flux is taken from SOLAR2000 (S2K) v.2.36 
model of Tobiska (2004) and scaled accordingly to the helio- 
centric distance of the comet at the time of observation. The 
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Figure 1: Photodissociative excitation cross sections for the 
production of 0( X D) from H2O and CO are taken from Hueb- 
ner et al. (1992). The photodissociation of CO2 for the pro- 
duction of 0( 1 D) are taken from Jain and Bhardwaj (2012). 
The photodissociation cross section of C0 2 producing 0( X S) 
is calculated using the yield suggested by Huestis et al. (2010) 
and total absorption cross section, ifc- represents the cross sec- 
tion value for the production of 0( X S) from H2O at 1216 A 
assuming 0.5% yield. For comparison the total photoabsorp- 
tion cross section of H 2 taken from Huebner et al. (1992) 
is also shown. The solar flux taken from SOLAR2000 (S2K) 
model on 9 March 1997 (solar minimum condition; solar radio 
flux F10.7 = 74.7 x 10~ 22 J s" 1 m" 2 Hz" 1 ) is shown with 
scale on the right side y-axis. 
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Figure 2: Calculated radial profiles for major production 
mechanisms of 0( 1 S) along with the total production profile 
for the abundances of 6% CO2 and 24% CO relative to H 2 
production rate of 8.3 x 10 30 s _1 . hi/ : solar photon and e t h 
: thermal electron 
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Figure 3: Calculated radial profiles for the photodissociation 
of CO2 producing 0( 1 S) at different wavelength bands for the 
abundances of 6% CO2 and 24% CO relative to H2O produc- 
tion rate of 8.3 x 10 30 s _1 . 



electron temperature profile required to calculate dissociative 
recombination rates is taken from Lovell et al. (2004). Bhard- 
waj and Raghuram (2012) have found that the yield of 0( 1 S) 
in the photodissociation of H 2 at solar H Ly-a can not be 
more than 1%. In the present study we have taken this yield 
value as 0.5%. The impact of this assumption was discussed 
in our previous work (Bhardwaj and Raghuram, 2012). The 
photodissociative excitation cross section for CO2 producing 
0( 1 D) is taken from Jain and Bhardwaj (2012). The photodis- 
sociative excitation cross sections for the production of 0( 1 D) 
and 0( 1 S) from H2O, CO2, and CO used in the model are pre- 
sented in Figure 1. The attenuation of solar radiation and so- 
lar UV-EUV generated photoclcctrons in the cometary coma 
are described in our previous works (Bhardwaj et al., 1990; 
Bhardwaj, 1999, 2003; Bhardwaj and Haider, 1999; Raghu- 
ram and Bhardwaj, 2012). 

3 Results 

3.1 Production and loss of 0( 1 S) 

The calculated 0( 1 S) volumetric production rate profiles for 
major production processes are presented in Figure 2. The 
photodissociation of CO2 is the major production process of 
0( 1 S). Above cometocentric distance of 1000 km, the pho- 
todissociative excitation of H2O is also an equally important 
production source of 0( 1 S). Photodissociative excitation of 
CO is the next significant production mechanism in produc- 
ing 0( 1 S). Since no cross section is reported in the literature 
for photodissociation of CO producing 0( 1 S), we have taken 
the photo-rate for this process from Huebner and Carpenter 
(1979) and assumed that the formation of 0( 1 S) is similar to 
0( 1 D). This assumption results in the calculated 0( X S) profile 
below 100 km similar to that of 0( 1 D). However, this assump- 
tion does not make any significant impact on the calculated 



green line intensity, since photodissociation of C0 2 and H2O 
can produce 0( 1 S) an order of magnitude higher than that of 
CO in the inner coma. Above 10 4 km, the contribution from 
dissociative recombination reactions of H20 + and CO + to the 
total 0( 1 S) production is significant. The photodissociation of 
OH is a minor source of 0( 1 S) below 10 5 km radial distance. 

The calculated 0( 1 S) volumetric production rate profiles for 
photodissociation of CO2 in the different wavelength bands are 
shown in Figure 3. The cross section for photodissociation of 
C0 2 in the wavelength band 955-1165 A is higher by a few 
orders of magnitude compared to that at other wavelength re- 
gions (cf. Fig. 1). Moreover, in this wavelength band, the yield 
of 0( 1 S) in photodissociation of C0 2 tends to unity (Slanger 
et al., 1977; Lawrence, 1972), while the total absorption cross 
section of H 2 has a strong dip (cf. Fig. 1). Thus, solar pho- 
tons in this wavelength band can dissociate C0 2 and produce 
0( 1 S) very efficiently. The photons in the wavelength bands 
1165-1375 and 745-955 A make a smaller (<10%) contribution 
to the total 0( 1 S) production. The contribution of 1216 A so- 
lar photons to the 0( X S) formation is two orders of magnitude 
low because of the small absorption cross section of CO2 (~8 
x 10~ 20 cm 2 ). 

The calculated volumetric destruction rate profiles of 0( X S) 
are presented in Figure 4. The collisional quenching of 0( 1 S) 
by H2O is the dominant loss process at cometocentric dis- 
tances shorter than 300 km. Above 1000 km the radiative 
decay via [OI] 5577 A line emission is the main loss process 
for the 0( 1 S) atom. The radiative decay via [OI] 2972 A emis- 
sion is a minor loss process of 0( 1 S). 

3.2 Production and loss of 0( 1 D) 

The calculated volumetric production rate profiles of 
metastable 0( 1 D) for different formation mechanisms are 
shown in Figure 5. Between 100 and ^2 x 10 km, most of 
the 0( 1 D) (>90%) is produced via photodissociation of H 2 0. 
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Figure 4: Calculated radial profiles for the major loss mecha- 
nisms of the 0( 1 S) atom for the abundances of 6% CO2 and 
24% CO relative to H 2 production rate of 8.3 x 10 30 s _1 . 



Figure 5: Calculated radial profiles for the major production 
mechanisms of 0( 1 D) along with the total 0( 1 D) production 
rate profile for the abundances of 6% CO2 and 24% CO relative 
to H 2 production rate of 8.3 x 10 30 s" 1 . hz^ : solar photon. 



However, below 100 km, the photodissociation of CO2 is also 
an important source of 0( 1 D). Between 200 and 2000 km, the 
radiative decay of 0( S) makes a minor contribution in the 
formation of 0( 1 D). Above 10 4 km, the photodissociation of 
OH plays a significant role in the formation of 0( 1 D). Even 
though the relative abundance of CO in Hale-Bopp is high 
(~25%), the photodissociation of CO is not a potential source 
mechanism of 0( 1 D). The calculated 0( 1 D) photodissociation 
rate profile for photodissociation of CO shows a double peak 
structure, which is explained later. 

The wavelength-dependent production rates of 0( 1 D) in the 
photodissociation of H 2 are shown in Figure 6. The most in- 
tense line of solar UV spectrum, H Ly-a at 1216 A, produces 
maximum 0( 1 D) around 1000 km, while solar photons in the 
wavelength regions 1165-1375 and 1375-1575 A are responsi- 
ble for producing maximum 0( X D) at shorter radial distances 
of 200 and 50 km, respectively. Since the total absorption cross 
section of H 2 in the 1165-1575 A wavelength region is small 
(cf. Fig. 1), these solar photons are able to penetrate deeper in 
the coma and mostly get attenuated at shorter cometocentric 
distances by dissociating H 2 0. The 0( 1 D) formation rate by 
solar photons at other wavelengths is smaller by more than an 
order of magnitude. 

Similarly, the production rate of 0( 1 D) due to photodis- 
sociation of C0 2 calculated at different wavelength bands is 
shown in Figure 7. At radial distances <100 km, solar pho- 
tons in 1375-1585 A wavelength region is the main source for 
0( 1 D) formation. This is because the absorption cross sec- 
tion of H 2 has a strong dip around 1400 A (cf. Fig. 1) and 
the average absorption cross section values of H 2 and C0 2 
are nearly same in this wavelength region. Thus, solar pho- 
tons in this wavelength band are able to reach the innermost 
coma and produce 0( 1 D) by dissociating C0 2 . Since the cross 
section for production of 0( 1 D) in photodissociation of C0 2 
peaks in the wavelength band 955-1165 A, the solar photons 
of this region leads the production of 0( 1 D) above 500 km. 

The production rates of 0( 1 D) via photodissociation of CO 



in different wavelength bands are presented in Figure 8. The 
total absorption cross section of H 2 is around two orders of 
magnitude smaller below 115 A than at other wavelengths, so 
these high energy photons can travel deeper into the cometary 
coma (even below 100 km) almost unattenuated. Since the CO 
molecule offers a cross section (average ~2 x 10~ 20 cm 2 ) to 
these photons it leads to the formation of 0( X D) and C( 1 D) 
via photodissociation closer to the cometary nucleus. Between 
100 and 500 km, the solar photons in the wavelength region 
115-325 A produce maximum 0( X D) atoms via photodissoci- 
ation of CO. The dissociative excitation cross section of CO 
is maximum in the wavelength region 535-955 A (cf. Fig. 1), 
which results in the peak production of 0( X D) via photodis- 
sociation of CO at 1000 km. More details on the attenuation 
of solar flux in high water production rate comets are given in 
Bhardwaj (2003). 

The model calculated volumetric loss rate profiles of 0( 1 D) 
are presented in Figure 9. This figure depicts that the predom- 
inant destruction channel of O( D) in the inner coma (below 
3000 km) of comet Hale-Bopp is quenching by H 2 0, which 
results in the formation of two OH molecules. Above radial 
distance of 10 4 km, the radiative decay leading to the red- 
doublet emissions is the major loss for 0( 1 D) atoms. Quench- 
ing by C0 2 and CO are minor loss processes, about one order 
of magnitude smaller and hence is not shown. 

The calculated density profiles of 0( 1 S), 0( 1 D), and 0( 3 P) 
in comet Hale-Bopp along with parent species considered in 
our model are shown in Figure 10. The density of 0( 1 S) peaks 
around 500 km, while the density profile of 0( : D) shows a 
broad peak between 2000 and 5000 km. The calculated num- 
ber density profiles of 0( 1 D) and 0( X S) without collisional 
quenching processes are also presented in this figure (with 
dashed lines). This calculation clearly shows that collisional 
quenching can significantly reduce the 0( X S) and 0( 1 D) den- 
sities in the inner coma. The formation of 0( 3 P) below 200 
km is due to collisions between OH molecules. 
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Figure 6: Calculated radial profiles for the photodissociation 
of H2O producing 0( 1 D) at different wavelength bands for the 
abundances of 6% CO2 and 24% CO relative to H2O produc- 
tion rate of 8.3 x 10 30 s _1 . 



3.3 Forbidden emissions of atomic oxygen: 
[OI] 5577, 2972, 6300, and 6364 A 

The emission rates of [OI] 5577, 2972, 6300, and 6364 A 
are calculated by multiplying Einstein transition probabilities 
(A 5577 = 1.26 s-\ A 2972 = 0.134 s"\ A 630 o = 6.44 x 10~ 3 s" 1 , 
and A 636 4 = 2.17 x I0~ 3 s~ 4 ) with the densities of 0( 1 S) and 
0( 1 D) (see Bhardwaj and Raghuram (2012) for calculation 
details). The intensity of these line emissions along the line 
of sight is calculated by integrating the emission rates. The 
model calculated brightness profiles as a function of projected 
distance for these forbidden emissions along with the [OI] 6300 
A observations of Morgenthaler et al. (2001) made on 2 and 
5 March 1997 using Hydra and WHAM instruments, respec- 
tively, are presented in Figure 11. To show the collisional 
quenching effect, we also presented the calculated forbidden 
emission line intensities (with dotted lines) in Figure 11, by 
considering only radiative decay as the loss process of 0( X S) 
and 0( 1 D). The [OI] 2972 A emission profile is shown by tak- 
ing branching ratio of 5577/2972 as 10 as suggested by Slanger 
et al. (2006). The NIST recommended value for this ratio is 
16 (Wiese et al., 1996). 

The calculated percentage contributions of various processes 
involved in the production of metastable 0( 1 S) and 0( 1 D) at 
different projected distances are presented in Table 2. For 
6% relative abundance of CO2, photodissociation of CO2 is 
the major source of 0( 1 S) production rather than photodis- 
sociation of H 2 (cf. Fig. 2). So we varied the CO2 relative 
abundance to study the change in the contribution of CO2 
to the 0( X S) and 0( X D) production. Calculations presented 
in Table 2 depict that, for a 6% relative abundance of CO2, 
below 10 4 km projected distances, around 25 to 30% of 0( 1 S) 
production is via photodissociation of H2O, while 40 to 60% 
production is through photodissociation of CO2 . Though the 
relative abundance of CO in comet Hale-Bopp is high (^25%), 
the photodissociation of CO could contribute a maximum of 
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Figure 7: Calculated radial profiles for the photodissociation 
of CO2 producing 0( X D) at different wavelength bands for the 
abundances of 6% CO2 and 24% CO relative to H2O produc- 
tion rate of 8.3 x 10 30 s _1 . 



10% to the 0( 1 S) production. The dissociative recombina- 
tion of H20 + and CO + together can contribute 10% to the 
production of 0( 1 S), whereas photodissociative excitation of 
OH is a minor (<5%) source. At 10 5 km projected distance, 
the photochemical reactions mentioned in Table 2 all together 
contributing 60% of 0( X S) and remaining is contributed by 
dissociative recombination of O-bearing ions. When the abun- 
dance of C0 2 is reduced to 3%, below 10 4 km projected dis- 
tance, photodissociation of H2O (35 to 40%) and CO2 (30 to 
50%) contribute almost equally to the production of 0( 1 S). 

The major production process of 0( 4 D) is the photodisso- 
ciation of H 2 0, whose contribution is 60 to 80% below 10 4 
km projected distance (cf. Table 2). Around 10 4 km the pho- 
todissociation of OH is also a significant production source of 
0( 4 D) and contributes around 20%; but, in the inner coma the 
contribution of this process is small (<10%). Radiative decay 
of 0( 4 S) and electron recombination of H 2 + contribute less 
than 10% each. At 10 5 km projected distance, most (75%) of 
0( 4 D) is produced by photodissociation of OH and remaining 
is contributed by other reactions. The change in the relative 
abundance of C0 2 by a factor of 2, from 6% to 3%, does not 
affect the relative contributions of various sources of 0( : D) 
below 10 4 km projected distance. 



Table 3: Calculated percentage contributions for the major 
production processes of green and red-doublet emissions in 
the total observed projected field of view (2.4 x 10 5 km) on 
comet C/1995 01 Hale-Bopp with varying relative abundance 
of C0 2 . 

C0 2 (%) W + H2O W + OH k + C0 2 e" + C0 2 + e- + H 2 + OfSHOfD) hi/ + TO 
6 23 (48)' 4(35) 41(3) 8(0.5) 7(5) (7) 7(1) 

3 32 (50) 6 (36) 30 (2) 5 (0.5) 10 (8) (7) 10 (1) 

1 42 (50) 8 (37) 13 (0,5) 2 (0.5) 13 (8) (5) 13 (1) 

''The values in parenthesis are the calculated percentage contribu- 
tions for red-doublet emission. 
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Table 2: Calculated percentage contributions for the major production processes of 0( 1 S) and 0( 1 D) in comet Hale-Bopp 
with varying relative abundance of CO2 for 0.5 % 0( 1 S) yield. 

CO2 Production processes of 0(*S) and 0( 1 D) at three comctocentric projected distances (km) (%) 

(%) hi/ + H 2 to + OH hi/ + C0 2 hy + C0 OCSHOCD) e + H 2 + 

if if if if if if if if if if if if W if if if 10 2 10 3 10 4 if 10' if if if~~ 

6 25 (77)' 31 (76) 24 (49) 6 (7) 0.5 (6) 1 (8) 4 (33) 14 (75) 58 (7) 50 (4) 40 (2) 23(1) 6 (0.5) 7 (0.5) 8 (1) 5 (0.5) (6) (8) (8) (5) 2 (2) 3 (2) 9 (8) 7 (5) 

3 33 (82) 42 (80) 33 (49) 8 (8) 1 (7) 1 (9) 5 (34) 18 (77) 47 (4) 36 (2) 26 (1) 15 (0.5) 8 (0.5) 9 (1) 10 (1) 7 (0.5) (5) (6) (6) (4) 2 (2) 3 (2) 11 (8) 11 (5) 

1 49 (85) 57 (82) 42 (51) 10 (8) 1 (7) 2 (9) 7 (35) 24 (78) 27 (1) 17 (0,5) 11 (0,5) 7 (0.5) 12 (1) 13 (1) 13 (1) 9 (0.5) (4) (5) (5) (3) 4 (2) 5 (3) 15 (7) 14 (5) 

^The values in parenthesis are for the 0( 1 D). 




Figure 8: Calculated radial profiles for the photodissociation 
of CO producing 0( 1 D) at different wavelength bands for the 
abundances of 6% CO2 and 24% CO relative to H 2 produc- 
tion rate of 8.3 x 10 30 s _1 . 



Figure 9: Calculated radial profiles for major loss mechanisms 
of the 0( 1 D) atom for the abundances of 6% CO2 and 24% 
CO relative to H 2 production rate of 8.3 x fO 30 s _1 . 



For a 4' circular aperture projected field of view (~2.4 x 
10 5 km) on comet Hale-Bopp, which is similar to the 50 mm 
Fabry-Perot spectrometer observations of Morgcnthalcr ct al. 
(2001), the calculated percentage contribution of major pro- 
duction processes for the green and red-doublet emissions, for 
different relative abundances of CO2, are presented in Table 3. 
These calculations clearly suggest that in a comet which has 
been observed over a large projected area, the photodissocia- 
tion of H 2 and OH mainly (~ 80%) controls the [OI] 6300 
A emission, while the radiative decay of 0( 1 S) contributes 
a maximum value of 10% to the total red-doublet intensity. 
With 6% relative abundance of CO2, the [OI] 5577 A line emis- 
sion observed in the coma is largely (~40%) contributed by 
photodissociation of CO2, and photodissociation of H2O is the 
next significant (*~25%) production process. The other pro- 
duction processes, like dissociative recombination of ions, pho- 
todissociation of CO, OH, etc, together contribute less than 
30% to the [OI] 5577 A intensity. When the C0 2 abundance 
is reduced to 3%, both photodissociation of H2O and CO2 
are contributing equally (~30%) to the green line emission in- 
tensity In all these cases, in spite of CO relative abundance 
being high (~25%) in comet Hale-Bopp, the photodissociation 



of CO could contribute a maximum value of 10%. 

3.4 Green to Red-doublet intensity ratio 

In comets, the parent species of these atomic oxygen emission 
lines are assessed using the ratio of intensity of the green line to 
the sum of intensities of the red-doublet, which can calculated 

as 

-^5577 _ T green a greenN green /3 green 

4300 + ^6364 7' J ^^a I . e( iA r red(/363O0+6364) 

where r is the lifetime of excited species in seconds (t[0( 1 D)] 
ps 110 s and t[0( 1 S)] sa 0.7 s), a is the yield of photodisso- 
ciation (Hucbncr ct al., 1992), f3 is the branching ratio (,$6300 
= 0.75, /3 63 64 = 0.25, /J5577 = 0.90, and /3 2g72 = 0.10 (Wiese 
et al, 1996; Slangcr et al, 2011; Festou and Fcldman, 1981) 
of the transition, and N is the column density of cometary 
species in cm" 2 . Customarily, the observed G/R ratio of 0.1 
has been used to confirm the parent species of these oxygen 
lines as H2O in comets (Cochran, 1984, 2008; Morrison et al., 
1997; Zhang et al., 2001; Cochran and Cochran, 2001; Furusho 
et al., 2006; Capria et al., 2005, 2008, 2010). However, since 
no experimental cross section or yield for the production of 
0( 1 S) from H2O is available in the literature, this ratio has 
been questioned by Hucstis and Slangcr (2006). In our pre- 
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Figure 10: Calculated number density profiles of 0( 1 S), 
0( 1 D), 0( 3 P), and OH, along with those of H 2 0, CO, and 
C0 2 . The calculations are done for the abundances of 6% 
C0 2 and 24% CO relative to H 2 production rate of 8.3 x 
10 30 s" 1 . The dashed lines of 0( 1 S) and 0( 1 D) are the cal- 
culated densities without accounting the collisional quenching 
processes. 



vious work (Bhardwaj and Raghuram, 2012), by fitting the 
observed green line emission intensity in comet C/1996 B2 
Hyakutake, we suggested that the yield for photo-production 
rate of 0( 1 S) from H2O at solar H Ly-a can not be more 
than 1%. Our previous model calculation also demonstrated 
that the determined G/R ratio depends on the projected area 
observed over the comet. 

We calculated the G/R ratio profiles on comet Hale-Bopp on 
26 March 1997 by varying CO2 relative abundance from 6 to 3 
to 1% which are presented in Figure 12. For comparison, the 
G/R ratio profile calculated on comet Hyakutake (Bhardwaj 
and Raghuram, 2012) is also plotted in Figure 12. In comet 
Hyakutake the G/R ratio is constant up to 100 km projected 
distance, while in the case of comet Hale-Bopp it is constant 
even up to 1000 km. The flatness of the G/R ratio depends 
on the quenching rate of metastable 0( 1 S) and 0( 1 D) by H 2 
which is a function of water production rate of the comet. 
Thus, in comets where H 2 production rate is still larger than 
that of Hale-Bopp, the G/R ratio would be constant up to 
projected distances larger than 10 3 km. 
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Figure 11: Calculated [OI] 6300, 5577, and 2972 A line bright- 
ness profiles as a function of the projected distances for the 
abundances of 6% C0 2 and 24% CO relative to H 2 produc- 
tion rate of 8.3 x 10 30 s _1 . For comparison the observed in- 
tensities of 6300 A line emission by Morgcnthalcr et al. (2001) 
using Hydra (stars) and WHAM (filled circles) instruments on 
1997 March 2 (r h = 1.05 AU and A = 1.46 AU) and 5 (r h 
= 1.03 AU and A = 1.42 AU), respectively, are also shown. 
During the observation the field of view of Hydra and WHAM 
instruments are 1° and 45', respectively (Morgenthaler et al., 
2001). Dotted lines are the calculated intensities when colli- 
sional quenching is not accounted. 



Hale-Bopp and Hyakutake, respectively. This figure depicts 
that in comet Hale-Bopp all the 0( 1 S) atoms produced above 
1000 km radial distance emit 5577 A (or 2972 A) photons, 
while for 0( 1 D) the radiative efficiency is unity above 10 4 km. 
Since the lifetime of 0( 4 D) is higher by two orders of magni- 
tude than that of 0( 1 S), most of the produced 0( 4 D) in the 
inner coma get quenched by other cometary species (mainly 
by H 2 0) without emitting photons at wavelengths 6300 and 
6364 A. But in case of comet Hyakutake the radiative effi- 
ciency of 0( 1 S) and 0( X D) is unity above 100 and 1000 km, 
respectively. This calculation shows that in comets most of 
the green and red-doublet emissions are produced above the 
collisional-dominated region where the radiative decay is the 
dominant loss process for 0( 1 S) and 0( 4 D) atoms. 



3.5 Radiative efficiencies of 0( 1 S) and 0( 1 D) 
atoms 

The number density of 0( 1 S) and 0( 4 D) in the cometary coma 
is controlled by various production and loss processes at that 
radial distance. To understand the region of maximum emis- 
sion of green and red-doublet lines in the coma we calculated 
the radiative efficiency profiles of 0( 1 S) and 0( : D) in comets 
Hale-Bopp and Hyakutake by calculating the ratio of emis- 
sion rate to total production rate of respective species. The 
calculated radiative efficiency profiles of 0( 1 S) and 0( 4 D) are 
presented in Figure 13 with solid and dotted line for comets 



3.6 Excess velocities of 0( 1 S) and 0( 1 D) 

Solar photons having energy more than the dissociation 
threshold of cometary species impart the additional energy 
to the kinetic motions of daughter products. The mean ex- 
cess energy released in the ith dissociation process at a radial 
distance r can be determined as 
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Figure 12: Calculated green to red-doublet intensity ratio along 
projected distances for different CO2 relative abundance [CO2] and 
with 0.5% yield for 0( 1 S) production in the photodissociation of 
H2O. Zhang et al. (2001) observed average green to red-doublet in- 
tensity ratio was 0.2 for the slit projected size of 522 x 1566 km over 
comet Hale-Bopp on 28 March 1997, which is shown with a hori- 
zontal line. The vertical dotted line represents 1566 km projected 
distance on the cometary coma. For comparison the calculated 
G/R ratio profile with 1% C0 2 and 0.5% yield in comet C/1996 
B2 Hyakutake is also shown. 



where A is the wavelength of solar photon, A t /j is the threshold 
wavelength for the dissociation process, h is Planck's constant, 
and c is the velocity of light. er(A) is the dissociation cross 
section of the cometary species at wavelength A. ^>(A, r) and 
t(A, r) are the solar flux and the optical depth of the medium 
for the photon of the wavelength A at a radial distance r, 
respectively. 



Our model calculated mean excess energy profiles for the 
photodissociation of H2O, CO2, and CO forming 0( 1 S) and 
0( 1 D) are presented in Figure 14 with solid and dotted lines 
for comets Hale-Bopp and Hyakutake, respectively. Above 
3000 km radial distance, the calculated excess energies in dif- 
ferent photodissociation processes in both comets show a con- 
stant profile, because the optical depth in this region for pho- 
tons of different wavelengths is very small. These values are 
in agreement with the calculations of Huebner et al. (1992). 
However, at shorter radial distances the neutral density is 
higher and hence the wavelength dependent photodissociation 
is significant which causes different excess energy values. 

In comet Hale-Bopp the calculated mean excess energy in 
photodissociation of H2O producing 0( 1 D) shows a highest 
value of 5.6 eV at the surface of the nucleus and decreases to 
a minimum value of 0.7 eV at 50 km. Above 50 km the mean 
excess energy increases and becomes constant (2.12 eV) above 
3000 km. This is because of the formation of 0( 1 D) via the 
photodissociation of H2O is associated with the photons of dif- 
ferent energies and it also varies with radial distance as shown 
in Figure 6. At a given radial distance the mean excess en- 
ergy released in the photodissociation process is determined by 




Cometocentric Distance (km) 

Figure 13: The calculated radiative efficiency profiles of 0( 1 S) 
and 0( 1 D) on comets Hale-Bopp (solid lines) and Hyakutake 
(dashed lines). Radiative efficiency is the ratio of emission 
rate to the production rate. 

the mean of energies of different solar photons involved. The 
threshold energy for production of 0( 1 D) by dissociating H2O 
is 7 eV. Very close to the cometary nucleus (<50 km), photons 
of wavelength smaller than 115 A and in the wavelength band 
1375-1575 A determines the formation of 0( : D) (cf. Fig. 6). 
At this distance, most of 0( 1 D) is produced by the photons 
of low energy (7-9 eV) in the wavelength band 1375-1575 A, 
and a small amount of 0( X D) is produced by very high energy 
(>100 eV) photons which results in the mean excess energy 
of about 2-5 eV. But around 50 km, the majority of 0( : D) 
production is determined by the photons of low energy 7 to 12 
eV (955-1575 A wavelength band) and the contribution from 
photons of wavelength below 115 A is very small. This causes 
the minimum value of mean excess energy 0.7 eV at this radial 
distance. 

Between 50 and 300 km, the increase in the excess energy 
is due to the production of 0( X D) atoms by photons of wave- 
length bands 115-325, 955-1575 A, and solar H Ly-a. Though 
high energy photons (115-325 A) are also involved in this re- 
gion, the intense solar photon flux at H Ly-a (1216 A) gov- 
erns the majority of O^D) production and subsequently deter- 
mines the mean excess energy. The solar H Ly-a photons can 
provide the maximum excess energy of 3 eV in the photodisso- 
ciation of H 2 0. Above 1000 km more than 90% of the 0( X D) 
production is controlled by photons at 1216 A wavelength and 
the remaining from other wavelength bands (cf. Fig. 6), which 
results a constant value of mean excess energy of 2.12 eV. 

Similarly, the mean excess energy calculated in the pho- 
todissociation of CO2 producing 0( X D) can be explained based 
on the wavelength dependent photon attenuated profiles pre- 
sented in Figure 7. The threshold energy for the 0( X D) pro- 
duction in photodissociation of CO2 is 7 eV and for 0( 1 S) it 
is 9 eV. At radial distances less than 100 km, the production 
of 0( 1 D) in photodissociation of CO2 is determined by the 
photons of low energy (average 8 eV) in the wavelength bands 
1375-1785 A and 955-1165 A, which results in low mean excess 
energy of ~1 eV. Above 100 km, photons of different energies 
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Figure 14: Calculated excess energy profiles of 0( 1 D) in pho- 
todissociation of H 2 0, CO, and CO2 and that of 0( 1 S) in 
photodissociation of CO2 on comets Hale-Bopp (solid lines) 
and Hyakutake (dashed lines). 

ranging from 7 to 16 eV (cf. Fig. 7) causes the mean excess 
energy of ~4 eV. The calculated mean excess energy profiles 
in the photodissociation of CO2 producing 0( 1 S) and 0( 1 D) 
are not similar. This is because the 0( S) production occurs 
via photodissociation of CO2 in the wavelength band of 800 
to 1300 A (photons of 10-15 eV), whereas 0( 1 D) can be pro- 
duced by photons of wavelength less than 800 A (>15 eV) (cf. 
Fig. 1). 

The threshold energy for the dissociation of CO produc- 
ing 0( : D) is 14.3 eV. Below 200 km the calculated maximum 
mean excess energy in the photodissociation of CO producing 
0( 1 D) is more than 100 eV. This is because the formation of 
0( 1 D) at these distances (cf. Fig. 8) is mainly determined by 
photons of wavelength less than 115 A (>110 eV) with some 
contribution from the wavelength band 115-325 A (40-110 
eV). Above 500 km, the formation of 0( 1 D) is mainly due to 
solar photons in the wavelength band 535-955 A (23-13 eV) 
which results in the maximum excess energy of 2.5 eV. 

4 Discussion 

The major difference between comets Hale-Bopp and Hyaku- 
take is the H2O production rate, which is larger by a fac- 
tor of 30 in the former. This difference in the H2O produc- 
tion rates result in a change in the photochemistry of 0( 1 S) 
and 0( 1 D) in the cometary coma. Due to the dense coma 
of comet Hale-Bopp, the attenuation of solar UV-EUV pho- 
tons on Hale-Bopp differs significantly from that in Hyaku- 
take. Moreover, the CO2 abundance in comet Hyakutake is 
smaller (<3% relative abundance) compared to that in Hale- 
Bopp (^6% relative abundance). The high H2O production 
rate in comet Hale-Bopp results in a larger collisional coma 
(radius few x 10 5 km) which is comparable to the scale length 
(~8 x 10 4 km) of H2O molecule. In the low production rate 
comets the collisional zone is smaller and photochemistry sig- 
nificantly differs. 



The photodissociation rates of H2O and CO2 for 0( 1 S) pro- 
duction differ by a factor of 20 (cf. Table 1). Hence, the 
major source of O^S) in the inner coma of comet Hale-Bopp 
is photodissociation of CO2 rather than photodissociation of 
H 2 0. Since the relative abundance of CO2 in comet Hyaku- 
take is 1%, the photodissociation of CO2 becomes an impor- 
tant source only near the surface of the nucleus (cf. Figure 
6 of Bhardwaj and Raghuram, 2012). The production peak 
of 0( X S) in comet Hyakutake is closer to the nucleus (<20 
km), whereas in comet Hale-Bopp it is between 100 and 1000 
km. Even when we reduced the CO2 abundance by 50% in 
Hale-Bopp, the peak production of 0( 1 S) in the inner coma 
is mainly controlled by photodissociation of CO2 and not by 
photodissociation of H2O. Hence, in a high water production 
rate comet a small relative abundance (~5%) of CO2, makes 
C0 2 as a potentially important source of 0( 1 S) compared to 
H 2 0. 

In comet Hyakutake, inside 10 5 km, the photodissociation 
of H2O is the major (more than 90%) production process of 
0( 1 D) formation and the contributions from other processes 
are very small. But in comet Hale-Bopp, since the H 2 pro- 
duction rate and CO2 relative abundance are higher, the solar 
photons of wavelength 955-1165 A, which are less attenuated 
by H2O, can travel deeper into the cometary coma and disso- 
ciate the CO2 to form 0( 1 D), which is not the case in comet 
Hyakutake. 

The radius of collisional coma, which is a function of to- 
tal gas production rate, in comets Hyakutake and Hale-Bopp 
differs by an order of magnitude. In comet Hyakutake quench- 
ing of 0( 1 S) by H 2 is the main destruction mechanism only 
close to the nucleus (<50 km) and radiative decay dominates 
at distances larger than 100 km. However, in comet Hale- 
Bopp collisional quenching is significant up to 500 km and 
only above that radiative decay is the major loss mechanism 
of 0( 1 S). Similarly, the collisional quenching radii of 0( X D) in 
comets Hyakutake (~10 3 km) and Hale-Bopp (~10 4 km) also 
differs by an order of magnitude. 

The 0( 4 D) density peak in comet Hale-Bopp is broader 
(2000 to 5000 km) than that in comet Hyakutake (200 to 
600 km). This change in the peak distribution of 0( : D) in 
the two comets is due to different H2O production rates and 
wavelength dependent photo-attenuation in the cometary co- 
mae. 

4.1 Comparison of model calculations with 
observations 

4.1.1 [OI] 6300 A emission 

Morgenthaler et al. (2001) observed [OI] 6300 A emission on 
comet Hale-Bopp on several days during February to April 
1997 using 4 different ground based instruments. Large aper- 
ture observations of 6300 A emission using WHAM and Hydra 
spectrometers are made for the field of view 1° and 45', which 
covers projected distances of 1.5 x 10 6 and 2.4 x 10 5 km on 
the comet, respectively. Our model calculated brightness pro- 
file of [OI] 6300 A emission shown in Figure 11 is consistent 
with these observations. The brightness profile of [OI] 5577 
A starts falling off beyond 1000 km, while the [OI] 6300 A 
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Table 4: The model calculated intensities of forbidden atomic 
oxygen emission lines on comet Hale-Bopp and the comparison 
of [01] 6300 A line with the observation of Morgenthaler et al. 
(2001) with 3% C0 2 and 24% CO. 



Date on 


r 


A 






Intensity (R) 




1997 


(AU) 


(AU) 


2972 A 


5577 A 


6364 A 


6300 A 














Calculated' 


Observed^ 


Mar 9 


0.999 


1.383 


34 


330 


1162 


3637 


2580-2922 


Mar 10 


0.992 


1.373 


36 


339 


1192 


3730 


2300-2649 


Apr 7 


0.920 


1.408 


15 


123 


1422 


4450 


2915-4964 


Apr 8 


0.923 


1.420 


13 


Ufi 


1400 


4379 


3057-3496 


Apr 9 


0.925 


1.431 


13 


111 


1380 


4323 


2920-3197 


Apr 10 


0.928 


1.444 


13 


103 


1358 


4248 


1579-1669 


Apr 13 


0.939 


1.484 


39 


372 


1271 


3296 


1451-1960 


Apr 14 


0.943 


1.497 


37 


361 


1240 


3878 


1575-2360 


Apr 16 


0.952 


1.526 


36 


339 


1179 


3688 


2335-2974 



*The calculated average surface 
jected distance of 2.5 x 10 5 km. 
§The upper and lower limits of 
Morgenthaler et al. (2001). 



brightness over the observed pro- 
[OI] 6300 A intensity observed by 



profile remains constant up to 5000 km. The flatness in the 
calculated surface brightness profiles mainly depends on the 
collisional quenching of metastable species which is a function 
of H2O production rate. The calculated green and red-doublet 
emission intensities (dotted lines in Fig. 11) when radiative de- 
cay is considered as the only loss mechanism shows the role of 
collisional quenching. Since the lifetime is relatively larger, the 
0( 1 D) is substantially quenched by H2O in the inner coma. 
Hence, below 1000 km, the calculated [OI] 6300 A emission 
intensities differ by a factor of 5. 

We also calculated the [OI] 6300 A emission intensity for 
a circular aperture of 4' diameter on different days of March 
and April 1997 similar to the observation conditions of Mor- 
genthaler et al. (2001). The calculated intensities of different 
atomic oxygen emissions are presented in Table 4 along with 
the [OI] 6300 A emission intensities observed by Morgenthaler 
et al. (2001). Our calculated intensities are higher by a factor 
of 1.5 to 2.5 compared to the observation and also vary for dif- 
ferent days due to change in solar flux, H2O production rate, 
and heliocentric and geocentric distances. The observed [OI] 
6300 A intensity on comet Hale-Bopp on 7 April 1997 is found 
to vary by a factor of 1.6 in a span of less than 20 minutes, 
which is difficult to explain with the variation of heliocentric 
dependent water production rate. Similarly, the observed in- 
tensity values on the other days of observation also show large 
variation. The brightness during April 10 to 14 is consistently 
lower than during April 7 to 9. The variation in the observed 
intensity might be associated with spectral extraction process 
because of non uniform sensitivity of Fabry-Perot spectrome- 
ter (Morgenthaler et al., 2001), rather than the intrinsic vari- 
ation in the comet. 

4.1.2 Green to red-doublet intensity ratio 

Zhang et al. (2001) observed comet Hale-Bopp on 26 March 
1997 using a rectangular slit (1.06" x 3.18") when the comet 
was at a geocentric distance of 1.32 AU and heliocentric dis- 
tance of 0.92 AU. For this observation, the projected field of 
view on the comet was 522 x 1566 km. Our calculated G/R 



ratio with 3% relative abundance of CO2 and 0.5% yield of 
0( 1 S) is 0.21, which is consistent with the observed G/R ra- 
tio range (0.18-0.22) of Zhang et al. (2001). The calculated 
average G/R ratio, for a 4' circular aperture field of view with 
3% relative abundance of C0 2 for the different days of obser- 
vation presented in Table 4, is around 0.1. This shows that in 
a high water production rate comet the observed G/R ratio 
over a large projected distances (~10 4 km) can be around 0.1 
(cf. Fig. 12). However, the calculated contributions of different 
production processes for 0( 1 S) suggest that photodissociation 
of CO2 is more important source rather than the photodisso- 
ciation of H2O. Hence, in comets with sufficient CO2 abun- 
dances (>5%), the green line emission is largely controlled by 
photodissociation of CO2 and the derived G/R ratio over large 
cometocentric distances could be around 0.1. 



Table 5: The model calculated green and red-doublet emission 
intensities and the derived 0( X D) and H2O production rates 
for different slit dimensions. The calculations are done with 
Q(H 2 0) = 8.3 x 10 30 s" 1 for the relative abundances of 6% 
C0 2 and 24% CO at r h = 1 AU and A = 1 AU using solar 
flux on 10 April 1997 (solar minimum period : solar radio flux 
F10.7 = 74.7 x 10~ 22 J s- 1 mT 2 Hz" 1 ). 

Average intensity (R) T Production rate 



Slit dimension 
(Projected distance in km) 
2" x 2" (725) 
5" x 5" (1.8 x 10 3 ) 
10" x 10" (3.6 x 10 3 ) 
30" x 30" (1.1 x 10 4 ) 
1' x 1' (2.2 x 10") 
4' x 4' (8.7 x 10 4 ) 
10' x 10' (2.1 x 10 5 ) 



[OI] 6300 A 
18895 [83188]* 
18909 [68301] 
19021 [52977] 
15668 [29341] 
11785 [18793] 
5005 [6767] 
2351 [3056] 



. 5577 A 
7584 [9245] 
6723 [7584] 
5369 [5825] 
2963 [3118] 
1846 [1924] 
605 [624] 
263 [271] 



"QPFl" 

1.3 x 10 26 

8.1 x 10 2C 

3.3 x 10 27 

2.5 x 10 28 

7.6 x 10 28 
5.0 x 10 29 

1.4 x 10 3 ° 



Q[H 2 0]' 
3.7 x 10 2(i 
23 x 10 26 

9.3 x 10 27 
7.1 x 10 2S 
2.1 x 10 29 

1.4 x 10 30 
3.9 x 10 30 



G/R" 

0.30 (0.08) n 
0.21 (0.08) 
0.21 (0.08) 
0.14 (0.08) 
0.10 (0.08) 
0.09 (0.07) 
0.08 (0.07) 



*The values in the square brackets are the calculated intensities 
without accounting for collisional quenching of 0( X S) and 0( 1 D). 
^Intensity is averaged over the projected field of view, 1 R = IjL 
Photons s _1 cm -2 sr _1 ; *The branching ratio for the production 
of 0( X D) in the photodissociation of OH is taken as 0.357 (see 
Morgenthaler et al., 2001), while for the photodissociation of H2O 
producing 0( 1 D) it is 0.064 (This work). The branching ratio (0.81) 
for the production of OH in photodissociation of H2O is taken from 
Huebner et al. (1992). ** Green to red-doublet emission intensity 
ratio determined over the projected field of view. ^The calculated 
G/R ratio without collisional quenching. 



To evaluate the role of slit dimension in determining the 
G /R ratio we calculated green and red line intensities for var- 
ious slit sizes by keeping H2O, CO and CO2 production rates 
as a constant. These calculations are presented in Table 5. By 
varying the slit dimension from 2" x 2" to 10' x 10' the cal- 
culated G/R ratio over the projected cometary coma changed 
from 0.3 to 0.08. This result clearly shows that the G/R ratio 
depends not only on the photo-chemistry in the coma but also 
on the projected area observed for the comet. The calculated 
G/R ratio is a constant value (0.08) throughout the cometary 
coma when collisional quenching is neglected in the model. By 
doubling the CO2 relative abundance in the coma, the G/R 
ratio increases by 30% whereas the collisional quenching of 
0( 1 D) and 0( 1 S) can change its value even by an order of 
magnitude. 
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Besides the dimension of the slit used for observation, the 
projected area observed on the comet depends on geocentric 
distance of the comet. Hence in a comet, where the collisional 
coma is resolvable in the observation, the derived G/R ra- 
tio depends on the projected area and also on the collisional 
quenching of 0( 1 S) and 0( 1 D) in the cometary coma. Thus, 
we conclude that the observed G/R ratio of 0.1 is not a defini- 
tive benchmark value to verify H2O or CO2/CO as the parent 
sources of atomic oxygen visible emissions in comets. 

4.1.3 Width of green and red-doublet emission lines 

Cochran (2008) has found that the width of green line is 
higher than either of the red-doublet lines in the spectra of 
8 comets. The wider green line implies the higher mean veloc- 
ity of metastable 0( 1 S), which could be associated with dif- 
ferent production processes. Besides collisions with different 
cometary species, the mean velocity of 0( 1 S) in the cometary 
coma is determined by various production processes, and/or 
could be due to the involvement of photons of various energies 
in dissociating O-bearing species (Cochran, 2008). 

The observed width of forbidden line emission depends 
on the velocity distribution of radiating metastable oxygen 
atoms. We found that the excess velocity released in photodis- 
sociation H 2 in the unity radiative efficiency region is 2.1 eV 
(cf. Figure 14). If we assume that most of this excess energy 
is transfercd to kinetic motion of atomic oxygen then the max- 
imum mean velocity that can be acquired by the 0( 1 D) atom 
would be 1.6 km s . This velocity is consistent with values 
of 0.5 to 1.8 km s^ 1 derived by Cochran (2008) in 8 comets. 
This supports the idea that most of the red-doublet emission in 
cometary coma is governed by the photodissociation of H2O. 
The excess energy profiles shown in Figure 14 suggest that 
the 0( 1 D) produced in photodissociation of CO and CO2 will 
have higher velocity than that produced in photodissociation 
of H 2 0. The excess energy released in the photodissociation 
of CO and CO2 in the unity radiative efficiency region is 2.5 
eV and 4.1 eV, which corresponds to 0( X D) excess velocity of 
^3.7 km s _1 and 4 km s _1 , respectively. However, our calcu- 
lations suggest that CO and C0 2 together can contribute to 
a maximum of 10% to the red-doublet emission. The contri- 
butions of CO and CO2 in the wings of red-doublet lines are 
probable. 

In the case of green line emission, since there is no experi- 
mentally determined cross section or yield for the photodisso- 
ciation of H2O producing 0( 1 S), it is difficult to determine the 
mean velocity acquired by an 0( 1 S) atom in the the photoly- 
sis of H 2 0. The maximum excess energy that can be released 
in photolysis of H2O producing 0( 1 S) at solar H Ly-a is 1.27 
eV. Again, if we assume all the excess energy is transferred as 
kinetic energy of atomic oxygen in S state then the maximum 
excess velocity of 0( 1 S) would be 1.3 km s _1 . But in the case 
of photodissociation of CO2, the excess energy is 2.5 eV, which 
corresponds to a maximum 0( X S) velocity of 4.3 km s _1 . The 
dissociative recombination of ions H2 + , CO^ , and CO + can 
contribute a maximum of 30% in the production of green line 
emission. But the excess energy released in these recombina- 
tion reactions is very small (Rosen et a!., 2000, 1998; Seiersen 
et al., 2003). By assuming that the maximum mean velocity 



that can be acquired by O^S) via the dissociative recombi- 
nation processes is about 1 km s _1 , we found that the mean 
velocity of 0( 1 S) from all production processes is ~2 km s _1 . 
This value is consistent with the derived velocity range of 1.9 
to 3.1 km s- 1 for 0( 1 S) in 8 comets by Cochran (2008). 

Before coming to a broad conclusion, we suggest that one 
has to calculate the exact mean excess velocities of 0( 1 S) and 
0( 1 D) over the observed cometary coma, by accounting for 
all collisional processes and the mean excess velocity profiles 
of various species. Due to non availability of photon cross 
sections for some of the photodissociation processes, and un- 
certainties involved in the excess energy calculations for disso- 
ciative recombination reactions, our model is limited in deter- 
mining the exact line widths of green and red-doublet emis- 
sions. However, based on our model calculations on comets 
Hale-Bopp and Hyakutake, we suggest that involvement of 
multiple sources in the formation O^S) could be a potential 
reason for the higher line width of green emission compared 
to that of red-doublet emission observed in several comets. 

4.2 Effect of model parameters on the calcu- 
lated intensities 

4.2.1 Expansion velocity of neutrals 

As we mentioned earlier in the Section 2, we have used the 
velocity profile from the work of Combi et al. (1999) for calcu- 
lating the number densities of parent species H 2 0, CO2, and 
CO. Combi et al. (1999) have shown that there is an accelera- 
tion of neutrals in the inner coma due to the photolytic heat- 
ing (Combi et al., 1999; Colom ct al., 1999; Biver et al., 1997; 
Combi, 2002) and other processes (Harris ct al., 2002). To 
evaluate the impact of this acceleration on our model results 
we carried out calculations by taking a constant gas expansion 
velocity profile with the values 0.7 and 2.2 km s . By using 
a constant velocity profile of 0.7 km s _1 in the coma, rather 
than a radially varying velocity of Combi et al. (1999), the 
calculated intensities of green and red-doublet emissions are 
increased by 30% and 25%, respectively, which are still higher 
than the observation. By changing the constant gas expansion 
velocity from 0.7 to 2.2 km s _1 , the calculated intensities of 
atomic oxygen emission lines are decreased by ~50%. However 
using the Combi et al. (1999) velocity profile, our calculated 
[OI] 6300 A emission intensities over 4' circular aperture field 
of view are closer to the observation (cf. Table 4). Hence, the 
velocity profile of neutral species is an important input in the 
model that should be accounted in calculating the intensities 
of these forbidden emissions. 

4.2.2 Relative abundances of neutral species 

The water production rate in comet Hale-Bopp has been de- 
rived using emissions of direct and daughter products of H2O 
by different observers (Weaver et al., 1997; Colom et al., 1999; 
Schleicher et al., 1997; Combi et al., 2000; Dello Russo et al., 
2000; Woods et al., 2000; Morgenthaler et al., 2001; Harris 
et al., 2002; Fink, 2009). During the observation period of 
these green and red-doublet emissions (r/j of the comet was 
around 0.9 AU), Dello Russo et al. (2000) measured the H 2 
production rates using infrared emissions of water molecules 



12 



for different days. In this period, Combi ct al. (2000) de- 
rived the H2O production rate in this comet using H Ly-a 
emission. The difference between these two derived produc- 
tion rates is less than 20%. These observations found that 
around 1 AU the water production rate in comet Hale-Bopp 
was about ~1 x 10 31 s _1 . Similarly, the derived water pro- 
duction rates of Fink (2009) on 1997 March 3 was 6.1 x 10 30 
s _1 which is smaller than the Combi et al. (2000) derived rate 
by a factor of 1.5. Using visible emission of atomic oxygen 
Morgenthalcr ct al. (2001) derived the H 2 production rates 
by applying standard branching ratios of OH and H 2 0. These 
derived H2O production rates are higher by factor of 3 to 6 
compared to values determined from other observations. To 
assess the impact of H2O production rate on the calculated 
green and red-doublet emissions we increased its value by a 
factor of 5. With increase in H 2 production rate the model 
calculated surface brightness of green and red-doublet emis- 
sions over 4' circular field of view is increased by a factor of 
3. 

As demonstrated earlier in this paper, the role of CO2 is very 
significant in determining the green line emission intensity and 
subsequently the G/R ratio. During the observation period of 
these forbidden emission lines the CO2 is not observed in this 
comet. To evaluate the impact of CO2 we varied its relative 
abundance from 3 to 6%. We found an increase (25%) in the 
calculated green line emission intensity over the 4' circular 
aperture field of view whereas it is small (<5%) for red-doublet 
emission intensity. 

Based on infrared observations made near perihelion on 
comet Hale-Bopp, DiSanti et al. (2001) suggested that 50% of 
CO abundance present in the cometary coma is contributed 
by distributed sources. Bockclcc-Morvan et al. (2010) investi- 
gated the extended distribution of CO by probing Hale-Bopp 
between ~800 to ^20,000 km region using CO rotational line 
emissions (viz, CO J(l-0) and CO J(2-l)). Based on the ob- 
servation and radiative transfer modelling studies, Bockelee- 
Morvan et al. (2010) rejected the idea of an extended distri- 
bution of CO in Hale-Bopp. Since the contribution of pho- 
todissociation of CO to formation of 0( X S) and 0( 1 D) is less 
than 10%, no significant variation in the calculated intensity of 
green and red-doublet emissions is found by reducing the CO 
relative abundance by half. Hence, the involvement of CO in 
these oxygen forbidden line emissions is almost insignificant. 

Though OH column densities are determined using 3080 
A surface brightness profile, there are large uncertainties in 
photo-cross sections of OH in producing 0( 1 D) and 0( 1 S) 
(Huebner et al, 1992; Morgenthaler et al., 2001). The calcu- 
lated photo-rates for the production of 0( 1 D) via photodisso- 
ciation of OH, using theoretical and experimental cross sec- 
tions differ by about an order of magnitude (Huebner et al., 
1992). Morgenthaler et al. (2001) studied the effect of these 
cross sections in deriving the H2O production rates using 6300 
A surface brightness profile and found that on using the the- 
oretical OH photodissociative branching ratios of 0( 1 D), the 
derived H2O production rates are higher by a factor of 3-6, 
than those determined based on experimental branching ra- 
tios of Nee and Lee (1984). The photodissociation of OH 
influences the calculated green and red-doublet emission in- 
tensities significantly above 10 4 km (cf. Figs. 2 and 5, and 



Table 2). By changing photorates determined by Nee and 
Lee (1984) experimental cross sections (which are used in the 
model) with the rates derived based on theoretically calcu- 
lated cross sections of van Dishocck and Dalgarno (1984), we 
found a 40% decrease in the calculated slit-averaged brightness 
over the 4' circular aperture field of view for both green and 
red-doublet emissions. But the calculated 0( 1 S) and 0( : D) 
production rates along the radial distances are decreased by 
an order of magnitude above 10 4 km. Since OH is the dom- 
inant O-bearing species in the outer coma, the cross sections 
can affect the calculated the surface brightness of [OI] 6300 A 
at larger projected distances (>10 5 km). To fit the observed 
[OI] 6300 A emission in the outer coma Glinski et al. (2004) 
found it necessary to increase theoretical determined OH to 
0( 1 D) photorate by a factor of around 3. 

The chemistry model developed by Glinski et al. (2004) sug- 
gested that the collisions of 0( 3 P) with OH leads to the for- 
mation of O2. These calculations also showed that the O2 den- 
sities can be as high as 1% of H2O. We evaluated the change 
in green and red-doublet emission intensities by incorporating 
O2 in the model by taking its density profiles from Glinski 
ct al. (2004). No significant change (<5%) is found in the 
green and red-doublet emission intensities by including O2 in 
the model. This is because the other O-bearing species are 
several orders of magnitude higher in the inner coma. 

4.2.3 Effect of slit dimension on the derived 0( X D) 
production rate 

As a case study, for a fixed H2O production rate and CO and 
CO2 relative abundances, we calculated [OI] 6300 A emission 
intensity over a projected field of view for different slit dimen- 
sions. We then derived the 0( X D) production rate based on 
the calculated average [OI] 6300 emission intensity over the 
projected field view. These calculations are presented in Ta- 
ble 5. Since our model calculations are limited up to the pro- 
jected distances of 10 5 km (which is discussed in Section 4.3) 
we present the calculated intensities of [OI] 6300 and 5577 
A emissions for the slit dimension up to 10' x 10'. Though 
0( 1 D) is substantially produced in the inner coma via pho- 
todissociation, the collisional quenching by cometary species 
results in a very few [OI] 6300 A photons. The role of quench- 
ing in determining the [OI] 6300 A flux can be understood 
from the calculated values presented in Table 5. A large aper- 
ture observation is required, which covers the entire [OI] 6300 
A emission region, to derive the H2O production rate. The 
calculations presented in Table 5 suggest that by using large 
aperture slit the derived water production rate is closer to 
the actual production rate of H 2 0. Hence, to derive the water 
production rate using [OI] 6300 A, the slit dimension which 
covers a projected distance more than the scale length of H2O 
should be used. 

4.3 Limitations and future scope of the 
model 

The density of the species produced in the inner coma (radial 
distances less than 10 5 km) is mainly controlled by photo- 
chemical reactions. Above these distances the transport of 
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species starts becoming significant in determining the num- 
ber density of the calculated species. Our model calculations 
are based on photochemical equilibrium condition and is for a 
collisional coma. Hence, model results presented at distances 
beyond 5 x 10 5 km are not as reliable as the values in the 
inner coma. Moreover, above these radial distances the chem- 
ical lifetimes of neutral species are significantly altered by the 
solar wind interaction through charge exchange and impact 
ionization processes. Also, we could not incorporate altitude 
distribution of dust density in our model calculations which 
can affect the calculated optical depth. Since our model is 
time independent and one dimensional it is difficult to explain 
the asymmetry in the observed [Of] 6300A emission intensity 
over the cometary coma. For determining the spectral width 
of green and red-doublet lines elaborated calculations are re- 
quired along with laboratory measured photodissociation cross 
sections. 

5 Summary and Conclusions 

We have recently developed a coupled chemistry-emission 
model for the forbidden visible emissions 5577 and 6300 A 
of atomic oxygen in comet C/1996 B2 Hyakutake (Bhardwaj 
and Raghuram, 2012). In the present paper we applied our 
model to a high (^30 times more than on Hyakutake) gas 
production rate comet C/1995 01 Hale-Bopp in which these 
prompt emissions are observed in 1997 by Morgenthaler et al. 
(2001) and Zhang et al. (2001). The main results of our model 
calculations on comet Hale-Bopp are summarized as follows. 

1. Below cometocentric distance of 10 3 km, photodissocia- 
tion of CO2 is the major production mechanism of 0( 1 S). 
Between 10 3 and 10 4 km, the contributions from the pho- 
todissociation of CO2 and H2O are nearly equal. Above 
2 x 10 4 km several other processes are also significant to 
the 0( 1 S) production. 

2. Mainly the solar photons in 955-1165 A wavelength band 
contribute to the production of 0( 4 S) in photodissocia- 
tion of CO2. This is because the yield of 0( 4 S) in CO2 
photodissociation reaches a maximum in this wavelength 
region. 

3. Since the cross section of photodissociation of CO2 for the 
production of 0( 4 S) is more than two orders of magnitude 
larger than that of H2O, even a small amount (few percent 
relative abundance) of CO2 can make it an important 
source of the 0( 4 S). 

4. Quenching by H 2 is the main loss mechanism for 0( X S) 
at radial distances below 300 km; above 10 3 km radiative 
decay via 5577 A emission is the dominant destruction 
mechanism. 

5. Inside 10 5 km, the main production mechanism of 0( 4 D) 
is photodissociation of H2O; but, in the innermost part 
of the coma (<100 km) the photodissociation of CO2 is 
also a significant source. 

6. For photodissociation of H2O, the peak 0( X D) production 
occurs via H Ly-a (1216 A), 1165-1375 A and 1375-1575 



A wavelength bands at cometocentric distances of 1000, 
200, and 50 km, respectively. Solar photons at all other 
wavelengths produce 0( 4 D) with one or more orders of 
magnitude smaller efficiency. 

7. Below 100 km, solar photons in the wavelength band 
1375-1585 A mainly produce 0( : D) by photodissociation 
of CO2. The contribution from other wavelength bands 
is significant above cometocentric distances of 200 km. 

8. The major destruction mechanism of 0( 4 D) up to 3000 
km cometocentric distance is quenching by H2O; above 
5000 km radiative decay takes over. 

9. In comet Hale-Bopp the 0( X D) density peaks occurs be- 
tween 10 3 and 10 4 km, while for 0( X S) the peak is around 
500-1000 km. 

10. The radiative efficiency of 0( 4 S) and 0( 4 D) atoms in 
comet Hale-Bopp are unity above 10 3 and 10 4 km, re- 
spectively. In comet Hyakutake these distances are 10 2 
and 10 3 km, respectively. 

11. The model calculated green to red-doublet emission in- 
tensity ratio is consistent with the observation of Zhang 
et al. (2001). 

12. Collisional quenching can change the G/R ratio by an 
order of magnitude, whereas doubling the relative abun- 
dance of CO2 increases its value by maximum of 30%. 

13. To accurately measure the H2O production rate in 
cometary coma, a slit dimension which covers a projected 
distance more than the scale length of H2O is preferred 
to cover the entire [OI] 6300 A emission region. 

14. The model calculated [OI] 6300 A emission intensity pro- 
file as a function of projected distance is in agreement 
with the observation of Morgenthaler et al. (2001). The 
model calculated surface brightness averaged over a 4' cir- 
cular aperture field of view is higher by a factor of 1.5 to 
2 compared to the observation. 

15. The calculated mean excess velocity of 0( X D) and 0( 4 S) 
atoms in the region of unity radiative efficiency is ~1.6 
and ^2 km s _1 , respectively, which is consistent with the 
range of velocities observed by Cochran (2008) in several 
comets. 

16. Based on our model calculations for comets Hyakutake 
and Hale-Bopp, we conclude that [OI] 6300 A emission is 
mainly controlled by the photodissociation of H 2 0, while 
the [OI] 5577 A emission line is contributed by both H2O 
and CO2. Since 0( 4 S) production is associated with dif- 
ferent molecules, whereas the O^D) production is mainly 
from H2O, the width of the green line will be higher than 
that of the red-doublet lines. 

With a high H2O production rate, comet Hale-Bopp pro- 
vided a large gaseous environment, which has not been seen in 
previous comets. Since the apparition was at small geocentric 
distances, the giant cometary coma has provided a laboratory 
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for investigating several collisional-driven effects. These col- 
lision driven processes are very important in determining the 
distribution of cometary excited species in the coma, which 
manifests into the emissions of the cometary coma. 
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